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Abstract 

At the end of 2010, the CERN Large Hadron ColUder started operation with heavy ion beams, col- 
liding lead nuclei at a centre-of-mass energy of 2.76 TeV/nucleon and opening a new era in ultra- 
relativistic heavy ion physics at energies exceeding previous accelerators by more than an order of 
magnitude. This review sunmiarizes the results from the first year of heavy ion physics at LHC ob- 
tained by the three experiments participating in the heavy ion program, ALICE, ATLAS, and CMS. 
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1 Introduction 

After almost two decades of design and construction, the LHC accelerator fP| started operation at the 
end of 2009 with proton-proton (pp) collisions at 0.9 and 2.36 TeV centre-of-mass energy, reaching its 
current top energy of ^/s = 7 TeV in 2010. Heavy ion collisions (A+A) are an integral part of the LHC 
physics program. The first four-week run with Pb beams took place in November 2010 at a centre- 
of-mass energy of 2.76 TeV per nucleon. About 10 pib^^ of integrated luminosity were delivered to 
each of the three LHC experiments which take part in the heavy ion program, the two general purpose 
detectors ATLAS and CMS [3] as well as the dedicated heavy ion detector ALICE [4,^. This review 
summarises the results from the first year of ion physics at the LHC, focusing on global event properties, 
high pt and heavy quark physics, and the production of hidden heavy flavour mesons (J/i/A, T). 

2 The Physics of Relativistic Heavy Ion Collisions 

Relativistic heavy ion collisions make it possible to study the properties of strongly interacting matter at 
energy densities far above those of normal nuclear matter. General arguments based on the property of 
asymptotic freedom suggest, and quantitative calculations of lattice quantum chromodynamics confirm 
ISl 121 18 |, that QCD matter undergoes a transition from a hadronic gas to a quark-gluon plasma at a 
temperature Tc 160 MeV, corresponding to an energy density of Ec ~ 0.5 GeV/fm^. At small net baryon 
density, the transition is a smooth cross-over spanning a temperature range of 20 — 30 MeV, which means 
that the precise value of the pseudo-critical temperature Tc depends on the observable used to locate it. 
The quark-gluon plasma phase is characterized by a much reduced condensate of light quarks, reflecting 
the approximate restoration of chiral symmetry, and by screening of the chromo-electric force between 
heavy quarks, implying the absence of quark confinement. At high net baryon density, the confined and 
deconfined phases are thought to be separated by a first-order phase transition, beginning with a critical 
point whose position is not accurately known. 

A large body of experimental and theoretical research conducted over the past decades at AGS, SPS, and 
RHIC has led to the following "standard model" of heavy ion collisions at collider energies |l9l|T0|. The 
energy deposited in the mid-rapidity kinematic range is controlled by the density of gluons contained 
in the colliding nuclei at moderately small values of the Bjorken scaling variable x fTTl. At top RHIC 
collision energy the relevant range is {x) k, 10^^, while at the current LHC energy the value is {x) 
10"^ In this range of x the gluon density obeys a highly nonlinear evolution equation that describes 
the saturation of the rapidly growing perturbative gluon distribution at low virtuality lfT2l . The virtuality 
scale at which this saturation occurs is called the saturation scale Qg. At present collider energies and for 
heavy nuclei 2^ ~ 1 .5 - 2 GeV^ at RHIC and 2^ ~ 3 _ 4 GeV^ at LHC [ 13 1. 

The liberation of the dense gluon sea by interactions during the collision of two nuclei is thought to 
result in the formation of a dense, non-thermal QCD plasma with highly occupied gauge field modes, 
often called a "glasma" llT4ll . The glasma thermalizes rapidly through nonlinear interactions including 
plasma instabilities |[T5l[T6l . During the longest phase of the heavy ion reaction the QCD matter forms a 
nearly thermal quark-gluon plasma, whose evolution can be described accurately by relativistic viscous 
hydrodynamics because of its very small kinematic shear viscosity. The temperature of the plasma at the 
time of thermalization is approximately IT^; it cools by mainly longitudinal expansion until it converts 
to a gas of hadron resonances when its temperature falls below Tc. At that temperature the chemical 
composition of the produced hadrons gets frozen, but the spectral distribution of the hadrons is still 
modified by final-state interactions, which are modelled by a Boltzmann transport equation. 

The very small dimensionless ratio of the shear viscosity T] to the entropy density s of the QCD matter 
produced in heavy ion collisions {Anr] /s <2.5 [17|) was one of the main discoveries of the experimental 
program at RHIC lITOl . Such a small ratio requires that the interactions among the constituent quanta 
are extremely strong and the cross sections among them approach the unitarity limit fTSl. This picture 
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is supported by QCD lattice simulations, which show that the QCD matter in the temperature range 
Tc <T < 2Tc is highly nonperturbative as witnessed, e. g., by the large value of the trace anomaly of the 
energy-momentum tensor, = e — 3P f/ll. 

The strongly coupled nature of the quark-gluon plasma created in the heavy ion reaction is confirmed by 
its ability to "quench" jets. At RHIC, the two main observables revealing this property are the suppression 
of single hadron yields at high transverse momentum pr relative to pp reactions and the additional 
suppression of back-to-back emission of high-;?7- hadrons. The ratios quantifying these suppressions are 
known as Raa and Iaa, respectively. Compelling arguments related the suppression effects to the ability 
of the quark-gluon plasma to degrade the kinetic energy of a hard parton (quark or gluon) traversing it, 
either by elastic collisions or collisions followed by gluon bremsstrahlung. The large values of the energy 
loss rate of hard partons deduced from the RHIC data [ 19 1 are consistent with the low shear viscosity 
deduced from the collective flow pattern of low-momentum hadrons ||20l . 

Before the LHC turned on, it was not clear whether this quantitative model of the dynamics of relativistic 
heavy ion collisions and the properties of hot QCD matter smoothly extrapolates from RHIC to the one 
order of magnitude higher LHC energies. Questions of particular interest were whether the increase in 
the initial temperature of the quark-gluon plasma would result in a much more viscous expansion of the 
created matter, and whether the gluon saturation effects in the colliding nuclei would be so strong that 
they significantly reduce, not only the expected density of the glasma, but also the yield of particles at 
moderately large pj. Because of the much higher beam energy and the large acceptance of the LHC 
detectors, heavy ion collisions at the LHC were also expected to provide access to new observables 
measuring the pattern of jet quenching and extend our knowledge of the interaction of heavy quarks and 
heavy quarkonia with hot QCD matter. An extensive compilation of theoretical predictions made before 
the commencement of the LHC heavy ion program can be found in ||2TI . 



3 Global event properties 

Global event properties describe the state and dynamical evolution of the bulk matter created in a heavy 
ion collision by measuring characteristics of the vast majority of particles which have momenta below a 
few GeV/c, referred to as "soft" particles. They include multiplicity distributions, which can be related 
to the initial energy density reached during the collision, yields and momentum spectra of identified 
particles, which are determined by the conditions at and shortly after hadronization, and correlations 
between particles which measure both size and lifetime of the dense matter state as well as some of its 
transport properties via collective flow phenomena. 



3.1 Multiplicity distributions 

The most basic quantity, and indeed the one measured within days of the first ion collision, is the number 
of charged particles produced per unit of (pseudo)rapidity, dA/^^h/dy (dA'^ch/dr]), in a central, "head-on" 
collision. When the LHC heavy ion program was conceived and detectors had to be designed more 
than a decade ago, predictions for dA'ch/dr] were highly uncertain, ranging from below 1000 to above 
4500, because experimental results had to be extrapolated from light ions at low energy (Sulphur beams 
at ^/s = 20 GeV) by orders of magnitude in both ion mass and beam energy ll22l |23 l . With results 
from RHIC, the uncertainties where substantially reduced, with most predictions concentrating in the 
range dN^h/dr\ = 1000 - 1700 EH. The value finally measured at LHC, dA^^h/dr] 1600 d, was 
on the high side of this range. From the measured multiplicity one can derive a rough estimate of the 
energy density with the help of a formula first proposed by Bjorken Ii25il relating the energy density to 
the transverse energy: 
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Fig. 1: a): Charged particle pseudorapidity density dNch/drj per colliding nucleon pair (O.SA^part) versus centre of 
mass energy for pp and A+A collisions, b): dA^ch/dJ] per colliding nucleon pair versus the number of participating 
nucleons together with model predictions for Pb+Pb at 2.67 TeV 



where To denotes the thermalization time, R is the nuclear radius, and Ej/N 1 GeV is the transverse 
energy per emitted particle. The value measured at the LHC implies that the initial energy density (at 
To = 1 fm/c) is about 15 GeV/fm^ ll26l . approximately a factor three higher than in Au+Au collisions at 
the top energy of RHIC E7ll28ll29ll30l[T0i . The corresponding initial temperature increases by at least 
30% to r 300 MeV, even with the conservative assumption that the formation time Tq, when thermal 
equilibrium is first established, remains the same as at RHIC. 

The high multiplicity at LHC, together with the large experimental acceptance of the detectors, allow for 
a very precise determination of the collision geometry (impact parameter and reaction plane orientation) 
in each event. Events are classified according to "centrality", where e.g. the 0-5% centrality bin contains 
the 5% of all hadronic interaction events with the largest final-state multiplicity and therefore the smallest 
impact parameters. Correspondingly, the 90-100% centrality bin comprises the 10% of events with 
the smallest final-state multiplicity and the largest impact parameters. The experimentally measured 
centrality can be related analytically or via a probabilistic Glauber Monte-Carlo generator to the number 
of participating nucleons Npa^t which are contained in the nuclear overlap volume for a given centrality 
class and therefore take part in the collision |31]. Likewise one can extract the number of nucleon- 
nucleon collisions A'coU between the participants. For Pb-i-Pb at LHC, the centrality or Np^n resolution 
varies from about 0.5% for central (A'part ~ 400) to about 5% for peripheral (A^part < 10) events Il32ll33l 

m. 

The charged particle multiplicity per participant pair ||32l[33l[34l, dNch/dri/{0.5 Npan), is shown in Fig.[T] 
together with lower energy data |[35l l36l l37l for central A-i-A collisions (typically 0-5% or 0-6% cen- 
trality). Note, however, that even for equal centrality selections, the average number of participants or 
collisions (A'pan, A'^coii) are different at the various energies because of different nuclear size (Pb or Au) 
and the rising pp cross-section. Particle production is no longer compatible with a logarithmic depen- 
dence with y^, as it was true for the data up to top RHIC energy |37|, but follows a power law ^ s^-^^. 
Also the pp data are well described by a power law, however with a less steep dependence on energy 

The centrality dependence of particle production is compared in Fig.[Tj5 with the one measured at RHIC, 
which is normalised to the LHC result at A'part^ 350 by scaling it with a factor 2.14. The results from 
the three LHC detectors Il32l [33l |341 are in excellent agreement with each other (within 1-2%) and have 
been averaged in this figure using the prescription of |35], assuming conservatively that the systematic 
errors in A'part are fully correlated between experiments. Comparison to the averaged and scaled 200 GeV 
Au-i-Au data (from ||35]| . updated using more recent STAR |[36l and PHOBOS 13711 data) shows a remark- 
able similarity in the shape of both distributions. For peripheral collisions, however, both distributions 
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extrapolate towards respective values measured in pp inelastic collisions (A'part= 2) at 200 GeV and 2.76 
TeV and therefore start to separate because of the different energy dependence seen for pp and A+A in 
Fig.[T}i. 

The fact that the shape of the normalised multiplicity distribution varies little with energy had already 
been noticed at RHIC 1371. It was still a surprise that it stays almost constant up to TeV energies, because 
hard processes, which scale with the number of binary collisions A'coih could be expected to contribute 
significantly to particle production at LHC and lead to a steeper centrality dependence, as predicted by 
the two component (soft + hard) Dual Parton Model DPMJET ll38i (dotted line in Fig. [T]^). However, a 
strong impact parameter dependent shadowing of the nuclear parton distribution function can limit this 
rise with centrality and is responsible for the flatter shape seen in the two component model HIJING [351, 
which agrees better with the data (full line). Saturation physics based on the "Colour Glass Condensate" 
(CGC) description [40], an example PTl is shown by the dashed line in Fig.[T]), naturally predicts such 
a strong nuclear modification, as well as the strong rise of particle production as a function of ^/s seen 
in Fig. [T^. In saturation models, the energy dependence is predicted to be a power law with roughly 
the correct exponent, whereas the centrality dependence changes only weakly (logarithmically) with the 
saturation scale and therefore the beam energy. 

3.2 Identified particle spectra 

Particle production (tt, K, p. A,..) is a non-perturbative process and cannot be calculated directly from 
first principles (QCD). In phenomenological QCD inspired event generators, particle spectra and ratios 
are adjusted to the data in elementary collisions (pp, e^e^) using a large number of parameters. In heavy 
ion reactions, however, inclusive particle ratios and spectra at low transverse momentum, which include 
the large majority of all produced hadrons even at LHC energies (about 95% of all particles are below 1.5 
GeV/c), are consistent with simple descriptions by statistical/thermal Il42ll43]| and hydrodynamical ll44l 
models, where particle ratios are determined during hadronisation at or close to the QGP phase boundary 
("chemical freeze-out", see below), whereas particle spectra reflect the conditions somewhat later in the 
collision, during "kinetic freeze-out". 

\npp collisions at high energy, the transverse momentum distribution can be separated into an exponen- 
tially decreasing low pj part (< few GeV/c), arising from soft, nonpertubative processes, and a powerlaw 
high pt tail from hard QCD scattering and fragmentation. A defining characteristic of heavy ion colli- 
sions is the appearance of ordered motion amongst the emitted hadrons in the soft part of the momentum 
spectrum ||45ll46| . It is called collective flow and implies, in contrast to random thermal motion, a strong 
correlation between position and momentum variables (nearby particles have similar velocities in both 
magnitude and direction). Flow arises in a strongly interacting medium in the presence of local pressure 
gradients. Different flow patterns are observed in heavy ion collisions and classified in terms of their 
azimuthal angle (p dependence with respect to the reaction plane, i. e. the plane spanned by the beam 
direction and the impact parameter. The isotropic (or angle averaged) component is called radial flow. 
In the framework of hydrodynamic models, the fluid properties (viscosity, equation-of-state, speed of 
sound, ..) together with boundary conditions both in the initial state (collision geometry, pressure gradi- 
ents, ..) and in the final state (freeze-out conditions) determine the pattern of collective motions and the 
resulting momentum spectra d^N/d/jj-dcp. 

Fig. [2^ shows the transverse momentum distributions of identified particles in central Pb-i-Pb collisions 
at the LHC |47i|48J. The spectral shapes differ significantly from both pp at LHC and Au-i-Au at RHIC, 
most dramatically for protons at low pj. The characteristic mass-dependent blue shift generated by the 
radial flow leads to a strong depletion at low pj - the proton spectrum is almost flat between 0.5 and 1 .5 
GeV/c - and to a harder spectrum at high pj, where the p/Ti ratio reaches a value of close to unity around 
3 GeV/c. 

The data are compared to hydrodynamic calculations, which are normalised to the data individually for 
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Fig. 2: a): Transverse momentum spectra of identified particles for central Pb+Pb collisions. The lines are 
scaled results from a boost invariant hydrodynamic model with (VISHNU, dashed lines) and without (VISH2+1, 
full lines) rescattering in the hadronic phase, b): Kinetic freeze-out temperatures Tfo and average radial flow 
velocities (j3) extracted from identified particle spectra (tt, K, p) at LHC (Pb+Pb) and RHIC (Au+Au) for different 
centralities. 



each particle species in order to better compare the shapes (absolute particle yields and ratios, which 
are an external input to the hydro models, are discussed below). The dashed line (VISH2+1) is a pure 
hydro calculation El |50l ED f or E and Q. (0-20% centrality); the full hne (VISHNU) for n, K, p (0- 
5% centrality) includes additional final-state rescattering, calculated with the URQMD transport code, 
which increases the radial flow velocity and improves the agreement with the data |52|. The hydro model 
describes the data fairly well up to intermediate momenta. Around 2 GeV/c, first pions and later protons 
start to deviate from the exponential hydro slope, indicating a progressive decoupling of high momentum 
particles from the thermalized and flowing bulk matter. 

The radial flow velocity can be estimated by fitting the spectra with a hydrodynamically inspired function 
called a "blastwave" fit [53]. The resulting kinetic freeze-out temperatures Tf„ and average radial flow 
velocities (j8) are shown in Fig.|2|5 for different centrality selections [47]. Quantitative results from such 
fits have to be taken with caution, as they in general depend on the particle types and momentum ranges 
included in the fit. The two fit parameters {Tfo, (j8)) are strongly correlated as indicated by the confidence 
contours in Fig. [3^ (each confidence contour corresponds to a different centrality class). However, a 
relative comparison to a similar fit to RHIC data [30] shows for the most central collisions that the 
average flow velocity increases significantly at LHC, reaching about 0.65c - the leading edge of the 
fireball expands collectively at essentially the speed of light - and that the kinetic freeze-out temperature 
drops below the one at RHIC. 

Fig. |3^ shows the ratio of A/K^ as a function of pj for central and peripheral nucleus-nucleus collisions 
at LHC and RHIC as well as at ^/s = 7 TeV [54]. Like the p/TT ratio mentioned above, the A/K^ ratio 
is strongly enhanced compared to pp, reaching a maximum of ?a 1 .5 around 3 GeV/c. Both the maximum 
value as well as its position in pT change only modestly compared to RHIC, where these unusual baryon 
to meson ratios had been observed previously and were referred to as "baryon anomaly" ll55l l56l l57]l . 
However, unlike at RHIC, the enhancement for central collisions of roughly a factor of three relative to 
the pp value persists out to at least 6 GeV/c. Qualitatively such an enhancement is expected from radial 
flow (the heavier baryons are pushed further out in pT than the lighter mesons), as shown by the full line, 
which represents the A/K^ extracted from the hydro model VISH2+1 already shown in Fig. |2^. Quan- 
titatively, however, the measured shape is very different indicating a progressive decoupling of hadrons 
from pure hydrodynamic flow above 2 GeV/c (see Fig. [2^). Global observables like momentum spec- 
tra, particle ratios and anisotropic flow seem to be governed by collective and thermal processes below 
Pt < 2 — 3 GeV/c, whereas hard processes and fragmentation dominate above > 6 — 10 GeV/c. In the 
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Fig. 3: a): The A/Kg ratio as a function of transverse momentum for two centrality selections in nuclear collisions 
(PbH-Pb at LHC, AuH-Au at RHIC) as well as pp collisions at 7 TeV. The full and dashed lines show the correspond- 
ing ratios from a hydro and a recombination model (scaled by 0.85), respectively, b) Particle ratios measured in 
central Pb-nPb (squares) and pp (triangles) collisions at the LHC. The full lines are the predictions from a thermal 
model. 



intermediate momentum region, various mechanisms have been proposed to describe the interplay be- 
tween soft and hard processes, including coalescence (recombination) of constituent quarks [58], strong 
colour fields with enhanced baryon transport [59 ], and colour transparency [60 1. The dashed line shows 
A/Kg calculated in a recombination model ||58l where radial flow and jet-quenching have been adjusted 
to the measured data at LHC; it describes the shape of the distribution above 3 GeV/c very well but is 
slightly off in normalisation (the model in the figure has been scaled by 0.85). A quantitative understand- 
ing at intermediate momenta remains a challenge and will require additional experimental constraints, in 
particular high precision data on identified particles out to at least 10 GeV/c. 

3.3 Identified particle yields 

Bulk particle production can be very successfully described in the framework of the thermal (statistical) 
hadronisation model Ii42ll43l . It assumes that particles are created in thermal (phase space) equilibrium 
governed by a scale parameter T, interpreted as a temperature. Production of a particle with mass m is 
suppressed by a Boltzmann factor e^'"^^ . Conservation laws introduce additional constraints, like the 
baryochemical potential Hb which accounts for baryon number conservation. An additional parameter 
Ys is introduced to describe the observation that in some collision systems particles containing strange 
quarks are suppressed compared to the grand canonical thermal expectation. The temperature parameter 
is found in all high energy collisions (pp, e^e^, A+A) to be about 160 - 170 MeV, while js increases 
from 0.5 — 0.6 in pp to 0.9 — 1 in A-i-A 16111621 . The disappearance of strangeness suppression in nuclear 
collisions ||63ll64| . usually called strangeness enhancement, was one of the first signals predicted for the 
QGP ||65]| . and the fact that the bulk of all particles are produced in heavy ion reactions with thermal ratios 
to very good approximation (typically < 10-20%) is considered to be an essential and well established 
fact. 

Fig.[3|5 shows a number of particle ratios from central Pb-i-Pb collisions, together with the (particle anti- 
particle averaged) values measured in pp at 1 TeV [47 , 66], as well as the prediction from a thermal model 
using the canonical conditions expected for nuclear collisions at the LHC ||67ll . Within experimental 
errors, particles and anti-particles are produced in Pb-i-Pb at midrapidity in equal numbers, consistent with 
the very small value of the baryochemical potential Hb of the thermal model. Strangeness is enhanced in 
Pb-i-Pb compared to pp, by a factor which increases with strangeness content from about L3 in KJn to up 
to > 3 in n/ 71 (note that the pp data in Fig. [3]) correspond to a higher energy than the Pb-i-Pb data). All 
strange particle ratios in Pb-i-Pb are very well described by the model, implying that they are produced in 
accordance with fully thermal ratios (i.e. js = !)■ However, in a complete surprise, protons were found 
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Fig. 4: a): Local freeze-out volume as measured by identical pion interferometry at LHC compared to central 
gold and lead collisions at lower energies, b): The system lifetime (decoupling time) T/ compared to results from 
lower energies. 



to be strikingly off - too low by a factor of about 1.5- and well outside the usual precision of the thermal 
model. The p/n ratio, which would have been expected to increase from its value measured in pp, stays 
essentially unchanged or even decreases slightly. 

Before concluding that something is wrong or missing in current implementations of the thermal model, 
the data, which are still preliminary, and various corrections (e.g. feed-down from weak decays) have to 
be thoroughly checked. Also the proton ratios at RHIC should be revisited, as there are indications of a 
smaller, but still significant, tension between model fits and data [68]. If the anomaly persists and cannot 
be described with extensions to the standard thermal model 16911701 . it may be due to inelastic processes 
which can change the particle composition after hadronisation. This mechanism is already thought to 
be responsible for the abundance of resonances with a large width and therefore a large hadronic cross 
section (like p and K*), but hard to describe quantitatively because of partially unknown cross sections 
and the difficulty to account for multi-particle initial states in hadronic transport codes. 



4 Particle correlations 

4.1 Identical particle (HBT) correlations 

The freeze-out volume (the size of the matter at the time when strong interactions cease) and the total 
lifetime of the created system (the time between collision and freeze-out) can be measured by identi- 
cal particle interferometry (also called Hanbury-Brown-Twiss or HBT correlations) {7T\. For identical 
bosons (fermions), quantum statistics leads to an enhancement (depletion) for particles emitted close-by 
in phase space. This modifies the two-particle correlation function, measured in energy and momen- 
tum variables, and can be related via a Fourier transformation to the space and time distribution of the 
emitting source, i.e. the space-time hyper surface of last rescattering. 

Results from HBT correlation measurements are shown in Fig. |4]for central collisions from very low 
energies up to LHC as a function of the charged particle density dNch/dt] [72|. The total freeze-out vol- 
ume is given as the product of a geometrical factor and the radii measured in three orthogonal directions 
(called /?iong5^side> and /?out)> whereas the lifetime was estimated from the pair-momentum dependence 
of Riong- The locally comoving freeze-out volume is directly proportional to the particle multiplicity 
(Fig. |4p) and therefore increases by a factor two compared to top RHIC energy to about 5000 fm^. The 
system lifetime is proportional to the cube root of the particle density (Fig. |4|)) and increases by about 
30% to 10 fm/c. Incidentally, the freeze-out volume roughly extrapolates at low beam energy to the 
volume of a Pb nucleus (« 800 fm^) and the lifetime vanishes at zero particle density. The evolution 
from RHIC to LHC of the individual radius parameters (/?iong7^side;^out) as well as their pair momen- 
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turn dependence is in satisfactory agreement with the predictions of hydrodynamical models lIlTl 1721 . 
The stronger transverse flow generated during the hydrodynamic phase at LHC energy reduces the im- 
portance of pre-equilibrium flow, in comparison with RHIC, and thus makes the agreement with the 
hydrodynamical predictions less dependent on extraneous assumptions 1731 . 

4.2 Anisotropic Flow 

4.2.1 General considerations 

The nuclear overlap zone in collisions with non-zero impact parameter is not azimuthally symmetric but 
has an almond shape whose deformation changes with centrality. Consequently, the pressure gradients 
between the centre of the overlap zone and its periphery in an average collision vary with azimuth, being 
strongest in the direction of the reaction plane angle ^rp, which coincides with the direction of the 
minor axis of the almond. The developing collective flow is proportional to the pressure gradient and 
therefore strongest towards the reaction plane, leading to an anisotropic distribution dN/d(p of particles. 
Anisotropic particle distributions were first suggested in 11741 as a signal of collective flow in ultra- 
relativistic heavy ion collisions and the flow pattern is usually quantified via a Fourier expansion |[75l : 

The Fourier (or flow) coefficients v„ depend on px and are given by Vn{pT) = (cos [n{(p — *P„)]), where 
the brackets denote an average over particles in a given pT bin and over events in a given centrality class. 
In the above equations, n is the order of the harmonic, (p is the azimuthal angle of the particle, and *P„ 
is the angle of the spatial plane of symmetry of harmonic n, the plane which maximises the expectation 
value of v„ in each event. Note that in general the symmetry angles *F„ can be independent from each 
other and point towards different directions in each individual event. 

Several experimental methods exist to measure the symmetry plane angles - using the phi-asymmetry 
generated by the flow itself to find the directions ^„ event-by-event - and the coefficients v„ (which 
are usually averaged over events), via two- or many-particle correlation measures l76l l45l l46ll . The vi 
coefficient is called directed flow . It is most prominent near beam rapidity but is not discussed further 
here (see fTTl ). Most attention was given in the past to V2, the elliptic flow, which is very strong in 
non-central collisions and directly linked to the almond shape overlap zone (^2 ~ ^rp)- Higher-order 
harmonics have usually been neglected because they were expected to be small for symmetry reasons 
(this assumption turned out to be surprisingly poor because of fluctuations, see below). 

The elliptic flow magnitude increases continuously with ^/s from SPS to RHIC Il45ll46l . At top RHIC 
energy, V2 reaches a value compatible with the one predicted by hydrodynamics for a "perfect fluid", i.e. 
a fluid without internal friction and vanishing shear viscosity ll27ll28ll29ll30llT0l . The shear viscosity, 
usually quoted in units oir\/s (shear viscosity T] over entropy density s), for a good relativistic quantum 
fluid is of order fi/kB- Using gauge gravity duality (conformal field theory in Anti de Sitter space, or 
AdS/CFT) as a proxy for strongly coupled QCD, it has been conjectured that a lower bound exists on 
77/5, i.e. A7ir]/s > 1 (/j = = 1) l78l . This value is reached when the 't Hooft coupling tends to infinity, 
and the mean free path approaches the quantum limit, the Compton wavelength. In such a "perfect" fluid, 
pressure gradients are not dissipated away because momentum transport perpendicular to the gradient is 
restricted by the ultra-short mean free path. 

Bounds on the 77 /s ratio can be derived from the data in two different ways. Either one compares the mo- 
mentum dependent elliptic flow parameter V2{pt) with results from viscous hydrodynamics calculations, 
or one extracts a value of by fitting the centrality (i. e. transverse size) dependence of the average pj 
integrated elliptic flow. Both approaches yield a T] ratio at RHIC that is very close to the lower bound, 
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Fig. 5: a): Elliptic flow coefficients V2 and V3 as a function of transverse momentum compared to a hydro 
model fSTl using two different values for the viscosity Tj /s b): Integrated elliptic flow coefficients as a function of 
collision centrality compared to a hydro calculation 18211 using two different models for the initial state geometry. 

at most 3-5 times 1/4?! 11791 [8OI . The most advanced model-data comparison using the second method 
gives an even tighter bound: Anrj/s < 2.5 [17]. 

Assuming no or only small changes in the transport properties of the quark-gluon plasma between RHIC 
and LHC, hydrodynamic models predicted that the (particle type averaged) elliptic flow coefficient V2, 
measured as a function of pr, should change very little, whereas the momentum integrated flow values 
should either stay approximately constant or rise by at most 30% 121]. This prediction was quickly 
confirmed [83]: the pr differential flow for charged particles is essentially unchanged at LHC whereas 
the pt integrated flow increases by some 30%, a direct consequence of the stronger radial flow discussed 
above. The matter created at LHC therefore still behaves like the (almost) perfect liquid discovered at 



In order to make progress in the precision determination of rj /s, some major obstacles have to be over- 
come. An example on the experimental side are non-flow correlations introduced e.g. by jets and reso- 
nance decays. On the theory side, uncertainties arise from model assumptions made about the initial state 
pressure gradients and their event-by-event fluctuations. It was in fact only recently realized Il84ll85ll that 
the statistical nature of individual nucleon-nucleon collisions can lead to highly irregular shapes of the 
reaction zone and the corresponding initial energy/pressure distributions. These shapes fluctuate from 
one event to the next, even at a fixed impact parameter. The irregular pressure gradients show no symme- 
try with respect to the reaction plane and therefore induce higher harmonic flow patterns. Fluctuations 
also tilt the symmetry angle of V2 away from the geometrical reaction plane {^'2 / ^rp), leading to 
event-by-event fluctuations of the elliptic flow direction and magnitude. 

4.2.2 Charged particle anisotropic flow 

The elliptic (vi) and triangular (V3) flow coefficients for mid-central Pb-Pb collisions are shown in Fig.|5^ 
as functions of pr- The agreement between the LHC experiments |[86l l87l [88l |89l [90l [9T| is again 
remarkable. Making use of the large rapidity acceptance of the LHC detectors, the coefficients v„ were 
measured by methods that minimize the influence of non-flow correlations, which are dominantly short- 
range in rapidity (2P = trigger di-hardon correlations, EP = Event Plane, SP = two-particle cumulant with 
eta-gap; see ||45| ). The elliptic flow coefficient rises approximately linearly with pj to a maximum of 
V2 ~ 0.23 around 3 GeV/c. The corresponding asymmetry in dA'^/d<p is very large indeed: almost three 
times as many 3 GeV/c particles aie emitted in-plane compared to out-of-plane (see eq|2]). The coefficient 
then decreases, at first rapidly then more gradually, but stays finite out to the highest pj measured. Also 
the triangular flow coefficient is very significant out to about 10 GeV/c and similar in shape, reaching 
about half the value of V2 at the maximum. 
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Fig. 6: a): Elliptic flow V2 as a function of pj for identified particles compared with a hydrodynamical model ll94l 
|95l . b): Elliptic flow for mesons and baryons, scaled by the respective number of valence quarks n^, versus scaled 
transverse kinetic energy. 



At high pt, where spectral shapes and relative particle yields indicate that hadrons are no longer in 
local thermal equilibrium, the nonzero V2 is thought to arise from differential parton energy loss rather 
than collective flow. The energy loss should depend on the path length inside the matter and therefore 
correlate naturally with the orientation of the elongated reaction zone. 

At lower pj, the data are compared to a hydrodynamic calculation lISTll using two different values of the 
shear viscosity(47rT]/5 = and 1). The effect of a nonvanishing viscosity is clearly visible: The initial 
pressure gradients are dissipated and the resulting collective flow is reduced. The influence of viscosity is 
seen to be much stronger on V3, as expected. The triangularity corresponds to shorter wavelength density 
variations and therefore larger local pressure gradients, which are more sensitive to viscosity. While both 
values of the viscosity are compatible with the V2 results, only A7ir]/s = 1 describes V3. A larger value of 
the viscosity would fall significantly below the V3 data, but could in principle be accommodated by using 
a different model for the initial state which generates a larger pressure gradient. 

This is illustrated in Fig.jS]?, which shows the pj integrated flow as a function of centrality f86l. The el- 
liptic flow rises strongly towards peripheral collisions, following the increasing elongation of the reaction 
region. Because of event-by-event fluctuations V2 stays non-zero even for the most central collisions. The 
triangular flow shows little centrality dependence, as it is driven solely by shape fluctuations. The hydro- 
dynamic model fST] shown in Fig. ^ employs two different, frequently used sets of initial conditions. 
The energy density distribution in the transverse plane is taken either from a Glauber calculation (solid 
line) combined with Anrj /s = I, or from a saturation model (dashed line) combined with Anrj/s = 2. The 
latter predicts larger pressure gradients and consequently needs a larger viscosity to describe the V2 data. 
However, the saturation model significantly underpredicts V3. Glauber initial conditions can describe 
both V2 and V3 better, in particular, if flow fluctuations are taken into account (see [82] for details). 

The measurement of higher harmonics thus seems to be able to overcome the degeneracy between shear 
viscosity and initial conditions and give a tighter limit of 47iri/s < 2 [S2]. Recent measurements of 
higher harmonic flow at RHIC have led to a similar conclusion 19211931 . However, none of the currently 
used initial state models can perfectly describe all experimental flow observations. With the on-going 
effort in both experiment and theory, the new high precision flow data from LHC promise significant 
further progress toward the goal of a precision measurement of the viscosity and other properties of the 
quark-gluon plasma. 

The interference between V2 and higher harmonic flow coefficients Il86ll96l 191119711981 . which are signif- 
icant up to « = 6, also gives a most natural explanation for some unusual long range structures observed 
in two-particle correlation at RHIC ll84l . traditionally called the "soft near-side ridge" and the "away- 
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side cone". The interpretation of these structures was controversial [99] until higher harmonic flow was 
firmly established at LHC and simultaneously at RHIC. 

4.2.3 Identified particle flow 

As for radial flow, a most stringent test of the collective flow interpretation of azimuthal anisotropics is 
the characteristic dependence on particle mass. The LHC data for vj for various particles (n, K, p, S, 
are shown in Fig. [6^ IiIOOliIOI iI . To understand the origin of the characteristic mass splitting, seen also in 
V3 mOOII . one needs to keep in mind that collective radial flow tends to equalise the velocities of particles, 
not their momenta, and therefore shifts heavy particles out to higher px than light ones II102I . The effect 
of a given azimuthal flow asymmetry thus manifests itself at higher momenta for particles with a larger 
mass. The hydrodynamical model [95^,94], which incorporates this effect, describes the data very well 
for all particle species up to intermediate pr- It also predicted the observation that the mass splitting is 
larger at LHC than at RHIC as a consequence of the increased radial flow. 

The mesons (n, K) deviate from the predicted flow above 1.5 GeV/c, whereas the baryons follow the 
hydrodynamical curves, within the still large experimental errors, out to about 3 GeV/c. The different 
behaviour of mesons and baryons, also seen in RHIC data [2T. ^lE^. '291 EHl [TOl, has been interpreted 
as a sign of quark recombination or coalescence 1 103, ,104. 105 , 106], already mentioned above as a 
possible reason for the large baryon-to-meson ratios. In this picture, hadrons at intermediate pj of a few 
GeV/c are formed by coalescence of constituent quarks from deconfined quark matter, with the three- 
quark baryons receiving a proportionally larger share of the collective quark flow than the two-quark 
mesons. Indeed, RHIC data indicate a universal scaling of flow HI 0711 for all particles of mass m when 
plotting V2/nq versus {mj —m)/nq), where Uq is the number of constituent quarks {uq = 2 or 3) and 

{piT — m) is the transverse kinetic energy {mj = ^ p\ + nP-). 

The same scaled distribution is shown in Fig. |6j) II 1001 IIOIII for LHC data. No scaling is apparent 
for pt < 1.5 GeV/c, where all particles follow more closely the mass dependent pattern predicted by 
hydrodynamics (see Fig. [6^). It therefore seems that the empirically observed low pj "quark scaling" 
at RHIC is a coincidence, only approximately valid for a particular value of the radial flow and the 
corresponding mass splitting. At higher pj, quark scaling between mesons and baryons may hold at 
LHC within the current experimental errors, but is definitely less convincing so far. Reduced errors and 
a measurement of flow for the (j) , a meson with a mass comparable to the proton, are required before a 
strong case for the coalescence picture can be made with LHC flow results. 

5 Hard processes at high pr and high mass 

The increased energy of heavy ion collisions at LHC relative to RHIC leads to much larger cross sections 
for hard processes, i.e. those involving high momentum or high mass scales. Energetic quarks or gluons 
can be observed as jets or single particles with pj reaching 100 GeV/c and beyond. Similarly, high 
Pt photons, charmonium and bottonium states (i.e. the J/i/a and T families), and even the weak vector 
bosons W and Z are copiously produced. The details of production and propagation of these high pj 
probes can be used to explore the mechanisms of parton energy loss and deconfinement in the medium. 

Some of the salient questions to be explored in this way are: How do colour charge and mass of a parton 
influence the energy loss? The perturbative QCD formalism of radiative energy loss [108] predicts that 
the energy loss of a parton should be proportional to the Casimir eigenvalue of its colour charge, i. e. that 
gluons should lose energy at roughly twice the rate of quarks, and that the rate of energy loss of heavy 
quarks should be reduced by the so-called dead-cone effect [109]. What is the relative contribution of 
radiative and elastic mechanisms of energy loss? How does the energy loss depend on the thickness 
of the medium? Strongly coupled gauge theories generally predict a stronger dependence on the path 
length in the medium than perturbative QCD llllOII . Where does the lost energy go? Does it remain 
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inside the jet cone and manifest itself as a modification of the fragmentation function, or does it get 
rapidly thermalized in the medium, making the jet look like a jet in vacuum, but with reduced energy? 

There are several experimental techniques available to explore these questions at the LHC. One is to 
measure single particle cross sections and compare them to the equivalent cross sections in pp colli- 
sions. The nuclear modification factor Raa of single particles, which has been extensively studied at 
RHIC Il27ll28ll29l[30l . allows for direct comparison with the RHIC results and their extension to higher 
momenta. The other approach is to reconstruct jets directly and compare jet yields in Pb+Pb with jet 
production in pp, as well as to measure jet-jet and jet-charged particle energy and angular correlations. 
The large kinematic range and the excellent calorimetry of the LHC detectors facilitate jet reconstruction, 
especially for jet energies in excess of 50 GeV. 

5.1 Single particle spectra 

The single particle production rates at RHIC have shown a large suppression of hadrons in nuclear 
collisions relative to pp, whereas particles that do not interact strongly, e. g. photons, are not modified. 
The LHC can significantly extend the accessible pj range and allow the measurement of additional 
particles, such as the Z and W. The suppression effects of a given particle are typically expressed in 
terms of the nuclear modification ratio: 

d^NAA/dpTdr] 

"^^P^^- {Tj^)d^o,,/dpjd^' (3) 

where Naa and o^n represent the particle yield in nucleus-nucleus collisions and the cross section in 
nucleon-nucleon collisions, respectively. The nuclear overlap function (7aa) is the ratio of the number of 
binary nucleon-nucleon collisions, {Ncow), calculated from the Glauber model, and the inelastic nucleon- 
nucleon cross section (cJj^^ = (64 it 5) mb at y/sm = 2.76 TeV). In the absence of nuclear effects the 
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factor Raa is unity by construction. As observed at RHIC in 200 GeV Au-i-Au collisions [i29l l30l . the 
yield of 5 — 10 GeV/c charged particles is suppressed in the most central events by more than a factor of 
five. Instead of Raa one can also approximate the centrality dependence by measuring Rqp, the ratio of 
central over peripheral events. 

While the charged particle Raa is the measurement with best statistical and systematic precision Ullll 
11121 1113] nT4 l. it is also interesting to measure the nuclear suppression factor for individual particle 
species to distinguish the exact mechanisms of energy loss. Measurements exist for identified K, Kg, 
A ImlinSl, isolated photons imi, Z,W UnmUllIISlini, ^-mesons una, jets im, Ifml 
I119[ll22 i| and T (VTT\. The suppression factors of prompt and non-prompt 7/ i/a mesons, the latter being 
produced from B-mesons and identified by their displaced decay vertex, were measured separately. In 
addition, ATLAS has presented results for the relative suppression of charged particles in central versus 



peripheral events, Rcp 111231 . 

A summary of Raa measurements for different particle species is shown in Fig. [7] for the most central 
events. The inclusive charged particle Raa follows, up to about 10 — 15 GeV/c, the characteristic shape 
discovered at RHIC (left panel, full circles). The pronounced maximum at a few GeV/c, which is some- 
times attributed to initial or final state interactions in nuclei ("Cronin effect"), is at very high energies 
more likely to be yet another manifestation of collective flow. It is qualitatively described by the dashed 
line, which shows the Raa obtained by dividing the inclusive charged particle distribution calculated by 
viscous hydrodynamics [50 | for central Pb-i-Pb by the experimentally measured pp spectrum. This inter- 
pretation is also supported by the fact that the apparent "suppression" factor is slightly larger for kaons 
and significantly larger for the A, as expected from flow. The peak region is followed by a steep decline 
and a minimum, around 5 — 7 GeV/c, where the suppression reaches a factor of about seven, very similar 
to but slightly larger than the one measured at RHIC. 



Heavy quarks, as shown by the Raa of prompt D mesons (open squares) and non-prompt J/i//^ (from 
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Fig. 7: Nuclear modification factor Raa as a function of pr for a variety of particle species together with theoret- 
ical predictions. Experimental error bars correspond to the total error (statistical and systematic errors added in 
quadrature), a) Low momentum region pj < 20 GeV; b) Entire momentum range measured at LHC. The curves 
show the results of various QCD-based models of parton energy loss 11241 11251 11261 [T27l 11281 . For details, see 
text. 



the decay of bottom quarks, closed diamond) in Fig. [7] are almost as strongly suppressed as inclusive 
charged particles. A similar conclusion can be drawn from the measurement of leptons from heavy 
flavour decays [ 1 15|. This seems contrary to the expectation that gluons, which are the dominant source 
of inclusive charged particles at LHC, should suffer twice as much energy loss as light quarks and that, in 
addition, the energy loss of heavy quarks should be even less than that of light quarks because of the mass 
dependence of radiation ("dead-cone" effect 111091 ). The strong suppression found for hadrons containing 
c- and Zj-quarks confirms observations made at RHIC and may indicate that the energy loss rate depends 
less strongly on the parton mass than expected for radiative energy loss. Reasons for this behaviour 
could be nonperturbatively large elastic energy loss in the strongly coupled quark-gluon plasma or heavy 
meson formation within the medium 111241 . More data and a quantitative comparison with models will be 
required to see how the small, with current statistics not very significant, difference between light hadron 
and heavy quark suppression can be accommodated by theory. 

Above pt ^ S GeV/c, the suppression becomes universal for all particle species (with the possible ex- 
ception of the non-prompt J/i/Aoriginating from B-meson decays shown in the left panel). With increasing 
Pt, Raa rises gradually towards a value of 0.5 (see right panel), a feature which was not readily apparent 
in the RHIC data. Isolated photons and the Z boson are not suppressed, within the currently still large 
statistical errors. This finding is consistent with the hypothesis that the suppression observed for hadrons 
is due to final-state interactions with the hot medium. 

The observed rise of Raa with pj allows a better discrimination between competing models of energy 
loss than the rather flat high pj dependence observed at RHIC. The rise can be understood as a decrease 
of the parton fractional energy loss with increasing pj, reflecting the weak energy dependence of pQCD 
radiative energy loss on parton energy. At RHIC this trend is compensated by the softening of the 
underlying parton spectrum, whereas at LHC the spectrum stays hard up to the highest measured pj 
which remains much farther away from the kinematic threshold than at RHIC. 

The observed trend is semi-quantitatively described by several models implementing the perturbative 
QCD (pQCD) formalism for energy loss lfT24l [T25l [T26l [T27l [T28l . The rate of induced gluon radiation 
in pQCD is governed by the rate of transverse momentum broadening, encoded in the jet quenching 
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parameter q, of the trajectory of a hard parton in the quark-gluon plasma flOSI. The value of q is 
proportional to the density of gluons in the medium and grows approximately as the entropy density of 
the plasma, or the charged particle multiplicity in the final state of the heavy ion collision. It should thus 
be slightly more than twice as large in Pb+Pb collisions at LHC than at the top RHIC energy. 

The band of curves labelled "GLV" II 12411 in Figure [7] shows results of a pQCD calculation using light- 
cone wave functions for the produced hadrons and assuming static scattering centres corresponding to a 
gluon density dNg/dy ^ 2,800 and a QCD couphng constant a, » 0.3. The curve labeled "YaJEM-D" 
111251 is a Monte-Carlo implementation of pQCD jet quenching that accounts for the limited virtuality 
evolution of the jet in a finite medium. The curves labelled "HT-W" II126I are based on the higher- 
twist formalism and correspond to values of ^oTq = 1.0 — 1.4 GeV^ at the time of thermalization Zq. 
The curve labelled "HT-M" [1271 is a different implementation of the higher-twist formalism with the 
parameter choice <?o = 1-3 GeV^/fm at Tq = 344 MeV. The "JEWEL" model |[T28l uses a Monte-Carlo 
implementation of an in-medium parton shower including LPM radiation suppression. 

The various calculations apply models of the nuclear reaction that are of widely varying detail. The 
GLV model uses a schematic description of the integrated medium density based on the boost-invariant 
Bjorken model of hydrodynamics. The YaJEM-D model combines ideal two-dimensional hydrodynam- 
ics with a saturation model for the initial energy density, but has not adjusted the parameters to the 
measured global properties of Pb-i-Pb collisions at LHC. The HT-W and HT-M calculations use three- 
dimensional ideal and two-dimensional viscous hydrodynamics, respectively, with initial conditions cho- 
sen to reproduce the charged particle multiplicity measured at LHC. The JEWEL model uses a variant 
of the schematic Bjorken hydrodynamical model. 

It is an important theoretical challenge to identify the origin of the differences between the five calcu- 
lations, in particular, the much more rapid rise of Raa in the YaJEM-D and HT-M models. Systematic 
studies of the differences between various implementations of radiative parton energy loss in pQCD 
at RHIC energies have revealed large sensitivities to poorly controlled aspects of lowest-order, coUinear 
gluon radiation. These differences can be further enhanced by the use of different models for the medium 
evolution. Nevertheless, the quantitative success of several models in correctly reproducing the rise and 
saturation of Raa with pj suggests that the energy loss of the leading parton in a jet shower may be 
described by perturbative QCD radiation in a strongly coupled medium. 

5.2 Jets 

Studying the modification of fully reconstructed jets is a particularly useful tool for probing the properties 
of the hot quark-gluon plasma. Jets are formed by fragmentation from high pj partons as they propagate 
through the produced matter. Measuring the energy of fully reconstructed jets allows one to distinguish 
between energy redistribution among the leading parton and the remainder of the jet and energy dissipa- 
tion out of the jet into the thermal medium. One of the promising channels are dijets, in particular their 
transverse energy balance and azimuthal angle correlation. The energy dissipation into the medium can 
be studied by measuring the asymmetry in pr of dijets in heavy ion collisions as a function of centrality 
and by comparing them to data from pp collisions. Similarly, the measurement of the relative azimuthal 
angle of the two jet axes gives information about the degree of scattering of partons as they traverse the 
medium. 

The measurement of the dijet asymmetry = {pn — pti) / {pt\ + Pti), where "1" and "2" refer to the 
leading and subleading jet, respectively, was performed by both ATLAS 1 129| and CMS [ 130|. Events 
containing at least two jets, with the leading (sub-leading) jet having pj of at least 120 (50) GeV/c for 
CMS and at least 100 (25) GeV/c for ATLAS, were selected for further study. 

The most striking observation by both experiments is the large centrality-dependent increase of the im- 
balance in the energy of the two jets, as measured in the calorimeters (See Fig.[8]l. While their energies 
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were very different, the two jets were observed to be very close to back-to-back in the azimuthal plane, 
implying little or no angular scattering of the partons during their traversal of the medium II 1301 as shown 
in Fig. [8] 1.131 il. 

The distribution of particle momenta inside jets normalized to the jet energy is the same, within experi- 
mental uncertainties, to that of jets produced in pp collisions as shown in Fig.|9][ 132|. This suggests that 
most of the additional energy radiated by the leading parton inside the medium gets absorbed by the mat- 
ter, and the fragments observed within the jet cone are produced outside of the medium. Several model 
calculations 111331 1134111351 [T36l . which combine elastic and inelastic parton energy loss with deflection 
of radiated gluons by the medium, have been able to reproduce the increased energy asymmetry of dijets 
in Pb-i-Pb. 
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Fig. 9: Ratio of Pb+Pb to jet fragmentation function as a function of jet asymmetry Ay for central Pb+Pb events. 
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The strong effect of jet quenching is confirmed by the ATLAS measurement of the jet Rcp il31l|. The 
comparison of jet production cross section between the central and peripheral events indicates that for 
the most central events the energy loss appears to reduce the cross section by about a factor of two (see 



Fig. 10 1. This shows that the energy lost by the leading parton is not simply redistributed within the jet 
cone but lost from the reconstructed jet. 
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Fig. 10: Jet central-to-peripheral event suppression factor Rqp for different centraUties. 



To find the "missing energy", the calorimetric measurement in CMS was complemented by a detailed 
study of low pt charged particles in the tracker and by using missing px techniques. The apparent 
missing energy in highly unbalanced dijet events was found among the low pr particles, predominantly 
with 0.5 < Pt <2 GeV/c, emitted outside of the sub-leading jet cone 1 130|. 

5.3 Quarkonium suppression 

Heavy quarkonia are important probes of the quark gluon plasma since they are produced early in the 
collision and their survival is affected by the surrounding medium II137II . The bound states of charm and 
bottom quarks are predicted to be suppressed in heavy ion collisions in comparison with pp, primarily 
as a consequence of deconfinment ("melting") in the QGP [138|. The magnitude of the suppression for 
different quarkonium states should depend on their binding energy, with strongly bound states such as 
the T showing less or no modification II139II . 

However, J/\j/ production, the classical deconfinement signal, has confounded expectations and interpre- 
tations ever since the first nuclear suppression was measured with Oxygen beams at the SPS 111401 ; now 
attributed to cold nuclear matter effects rather than deconfinement 1 141 1. The "anomalous" suppression 
seen later with heavier beams [ 142 | turned out, surprisingly, to be rather similar in magnitude at SPS and 
RHIC. This could indicate suppression of only the high mass charmonium states y/' and Xc, which popu- 
late about 40% of the observed J/y/, and which should dissociate very close to or even below the critical 
transition temperature. Alternatively, it has been suggested that the increasing (with energy) J/i/A suppres- 
sion is more or less balanced by enhanced production via recombination of two independently produced 
charm quarks II1431 11441 . A resolution to this puzzle may come from measuring J/y/ production at the 
LHC where coalescence effects should be stronger because of the abundant charm production, and by 
comparing the suppression patterns of the J/xj/ and T families. 

A compilation of first LHC results on quarkonia production II121I I122]| for both J/y (a) and T (b) is 
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Fig. 11: Nuclear modification factor Raa as a function of centrality for J/y/ (right) and Y (left). 



shown in Fig. 1 1 as a function of centrality (Npart), together with data from RHIC II1451I146TI1471 . 



While errors are still large, and the overall amount of suppression at LHC remains qualitatively similar 
to RHIC, the detailed pattern is quite different and intriguing. The pr integrated Raa measured for 
the J/\j/ at forward rapidity (closed circles) of about 0.5 depends very little on centrahty and is almost 
a factor of two larger than the one measured at RHIC in central collisions, also at forward rapidity 
(open circles); the difference is smaller but still significant when comparing with RHIC midrapidity 
data (open squares). On the contrary, the high pj data at LHC (full squares), which is compatible with 
an independent Rcp measurement 111 191 . show a stronger suppression than the high pr RHIC results 
(open stars). While such a pattern would be unexpected in a pure suppression scenario, it is qualitatively 
consistent with the recombination model, which predicts substantial regeneration effects only at low 
transverse momentum. The T suppression (right panel) is very similar at RHIC and LHC. As only about 
50% of the observed T(1S) are directly produced, and the T(2S/3S) states seem to be more suppressed 
than the ground state f 148], the measured Raa is compatible at both RHIC and LHC with suppression 
of the high mass bottonium states only. 

However, there are still a number of unknown aspects that must be clarified before firm conclusions can 
be drawn. Foremost are cold nuclear matter effects, in particular shadowing or saturation of the nuclear 
parton distribution functions, which for gluons are very poorly constrained at the small x values relevant 
at LHC. Likewise, very little is known about final state absorption of quarkonia in nuclei, which may 
(or may not) be important at LHC. For this reason a p+Pb run at LHC, which will address these nuclear 
effects, is mandatory and anticipated for 2012. 



6 Conclusions 

The LHC has entered the field of ultra-relativistic heavy ion physics with an impressive first year per- 
formance, having benefited from many years of preparation, a strong and very complementary set of 
state-of-the-art detectors, and, crucially, decades of experience and progress made at previous machines. 
Global event properties were found to evolve smoothly from RHIC to LHC, to be qualitatively similar 
but quantitatively different, and to reflect the expected increases in energy density, volume, and lifetime 
of the hot and dense matter created in the collisions. The data have generally confirmed that also at LHC 
energy the bulk features imply the formation of a very strongly interacting, nearly perfect quark-gluon 
plasma hquid. The prominent presence of various collective flow patterns has opened a path forward to 
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precision measurements of its material properties, including shear viscosity, equation-of-state, and sound 
velocity. 

Concerning bulk properties, and with the possible exception of thermal particle production, where the 
unexpectedly low proton yield needs confirmation and clarification, the standard reaction model estab- 
lished at RHIC has proven itself to be predictive and detailed enough to permit quantitative extrapolation 
to much higher colUsion energies. 

For observables where the energy reach of LHC is unique, i. e. low-x parton physics and hard processes, 
already the first low luminosity run has revealed a wealth of significant, and at times surprising results. 
The energy loss of high momentum partons traveling through the matter was measured with several dif- 
ferent particle species as well as with fully reconstructed jets. While its magnitude is roughly in the 
expected range, the more differential analysis of parton type and momentum dependence, dijet corre- 
lations, and energy flow around the jet axis has revealed features which were not anticipated and may 
spur substantial refinements of our understanding and modelling of parton interactions with the hot QCD 
medium. 

The results on quarkonium production, where precision measurements of the T resonances are possible at 
LHC, in addition to those of the charmonium states, have yet to provide decisive information on the sup- 
pression mechanism. However, intriguing hints have been found in the charmonium and bottomonium 
suppression pattern, suggesting significant differences to the observations at RHIC. 

The near future of heavy ion running at LHC will bring a significant increase in integrated luminosity - 
akeady the second heavy ion run in November 201 1 has increased the statistics for hard probes by more 
than an order of magnitude - as well as new signals like y-jet and Z-jet correlations. The p+Ph run 
planned for 2012 should not only deliver the much needed comparison data for the heavy ion program 
but also far extend the kinematic reach for the study of initial-state effects in the nuclear parton wave 
function, such as saturation and other low-;ic physics. And another factor of two increase in energy is in 
sight when LHC will reach the design value after 2014. 

The combination of the on-going heavy ion program at RHIC, which includes high-luminosity mea- 
surements with upgraded detectors, the new higher energy data from LHC, and the demonstrated solid 
foundation in theory promises rapid progress in a quantitative understanding of the properties of the 
quark-gluon plasma. 
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